Finite temperature dynamics in a polarized sub-Ohmic heat bath: A hierarchical equations of motion-tensor train study.

The dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperatures is investigated by employing both analytical tools and the numerically accurate hierarchical equations of motion-tensor train method. By analyzing the features of nonequilibrium dynamics, we discovered a bifurcation phenomenon, which separates two regimes of the dynamics. It is found that before the bifurcation time, increasing temperature slows down the population dynamics, while the opposite effect occurs after the bifurcation time. The dynamics is highly sensitive to both initial preparation of the bath and thermal effects.

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy.

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Finite temperature dynamics of the Holstein-Tavis-Cummings model.

By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra.

Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study.

This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.

Finite-Temperature Hole-Magnon Dynamics in an Antiferromagnet.

Employing the numerically accurate multiple Davydov Ansatz in combination with the thermo-field dynamics approach, we delve into the interplay of the finite-temperature dynamics of holes and magnons in an antiferromagnet, which allows for scrutinizing previous predictions from the self-consistent Born approximation while offering, for the first time, accurate finite-temperature computation of detailed magnon dynamics as a response and a facilitator to the hole motion. The study also uncovers a pronounced temperature dependence of the magnon and hole populations, pointing to the feasibility of potential thermal manipulation and control of hole dynamics. Our methodology can be applied not only to the calculation of steady-state angular-resolved photoemission spectra but also to the simulation of femtosecond terahertz pump-probe and other nonlinear signals for the characterization of antiferromagnetic materials.

Dynamics of a Magnetic Polaron in an Antiferromagnet.

The t-J model remains an indispensable construct in high-temperature superconductivity research, bridging the gap between charge dynamics and spin interactions within antiferromagnetic matrices. This study employs the multiple Davydov Ansatz method with thermo-field dynamics to dissect the zero-temperature and finite-temperature behaviors. We uncover the nuanced dependence of hole and spin deviation dynamics on the spin-spin coupling parameter J, revealing a thermally-activated landscape where hole mobilities and spin deviations exhibit a distinct temperature-dependent relationship. This numerically accurate thermal perspective augments our understanding of charge and spin dynamics in an antiferromagnet.

Contracted description of driven degenerate multilevel quantum systems.

We formulate a contraction theorem that maps quantum dynamics of a multilevel degenerate system (DS) driven by a time-dependent external field to the dynamics of the corresponding contracted non-degenerate system (CNS) of lower dimension, provided transitions between each pair of degenerate levels in the DS have identical transition dipole moments. The theorem is valid for an external field of any strength and shape, with and without rotating wave approximation in the system-field interaction. It establishes explicit relations between DS and CNS observables, significantly simplifies numerical calculations, and clarifies physical origins of the field-induced DS dynamics.

Thermo-Field Dynamics Approach to Photo-induced Electronic Transitions Driven by Incoherent Thermal Radiation.

The effects of thermal light-matter interaction on the dynamics of photo-induced electronic transitions in molecules are investigated using a novel first principles approach based on the thermo-field dynamics description of both the molecular vibrational modes and of the radiation field. The developed approach permits numerically accurate simulations of quantum dynamics of electronic/excitonic systems coupled to nuclear and photonic baths kept at different temperatures. The baths can be described by arbitrary spectral densities and can have any system-bath coupling strengths. In agreement with the results obtained previously by less rigorous methods, we show that the excitation process obtained by the continuous interaction with the suddenly turned-on thermal radiation field creates a mixed ensemble having a nonnegligible component consisting of a superposition of vibronic eigenstates which can sustain coherent oscillations for relatively long times. The results become especially relevant for the dynamics of electronic transitions upon sunlight excitation. Analytical results based on time-dependent perturbation theory support the numerical simulations and provide a simple interpretation of the time evolution of quantum observables.

Two-dimensional fluorescence excitation spectroscopy: A novel technique for monitoring excited-state photophysics of molecular species with high time and frequency resolution.

We propose a novel UV/Vis femtosecond spectroscopic technique, two-dimensional fluorescence-excitation (2D-FLEX) spectroscopy, which combines spectral resolution during the excitation process with exclusive monitoring of the excited-state system dynamics at high time and frequency resolution. We discuss the experimental feasibility and realizability of 2D-FLEX, develop the necessary theoretical framework, and demonstrate the high information content of this technique by simulating the 2D-FLEX spectra of a model four-level system and the Fenna-Matthews-Olson antenna complex. We show that the evolution of 2D-FLEX spectra with population time directly monitors energy transfer dynamics and can thus yield direct qualitative insight into the investigated system. This makes 2D-FLEX a highly efficient instrument for real-time monitoring of photophysical processes in polyatomic molecules and molecular aggregates.

Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation.

Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore.

By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.

Dynamics of dissipative Landau-Zener transitions in an anisotropic three-level system.

We investigate the dynamics of Landau-Zener (LZ) transitions in an anisotropic, dissipative three-level LZ model (3-LZM) using the numerically accurate multiple Davydov D2Ansatz in the framework of the time-dependent variational principle. It is demonstrated that a non-monotonic relationship exists between the Landau-Zener transition probability and the phonon coupling strength when the 3-LZM is driven by a linear external field. Under the influence of a periodic driving field, phonon coupling may induce peaks in contour plots of the transition probability when the magnitude of the system anisotropy matches the phonon frequency. The 3-LZM coupled to a super-Ohmic phonon bath and driven by a periodic external field exhibits periodic population dynamics in which the period and amplitude of the oscillations decrease with the bath coupling strength.

Dynamics of the spin-boson model: The effect of bath initial conditions.

The dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac-Frenkel time-dependent variational approach with the multiple Davydov D1 Ansatz in the interaction picture. The validity of our approach is carefully checked by comparing the results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for the spectral exponent s < sc, there exists a transition from coherent to quasicoherent dynamics with increasing coupling strengths. For sc < s ≤ 1, the coherent to incoherent crossover occurs at a certain coupling strength and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has a considerable influence on the dynamics.

Exciton Dynamics and Time-Resolved Fluorescence in Nanocavity-Integrated Monolayers of Transition-Metal Dichalcogenides.

We have developed an ab initio-based, fully quantum, numerically accurate methodology for the simulation of the exciton dynamics and time- and frequency-resolved fluorescence spectra of the cavity-controlled two-dimensional materials at finite temperatures and applied this methodology to the single-layer WSe2 system. Specifically, the multiple Davydov D2 Ansatz has been employed in combination with the method of thermofield dynamics for the finite-temperature extension of accurate time-dependent variation. This allowed us to establish dynamical and spectroscopic signatures of the polaronic and polaritonic effects as well as uncover their characteristic time scales in the relevant range of temperatures. Our study reveals the pivotal role of multidimensional conical intersections in controlling the many-body dynamics of highly intertwined excitonic, phononic, and photonic modes.

Plenty of Room on the Top: Pathways and Spectroscopic Signatures of Singlet Fission from Upper Singlet States.

We investigate dynamic signatures of the singlet fission (SF) process triggered by the excitation of a molecular system to an upper singlet state SN (N > 1) and develop a computational methodology for the simulation of nonlinear spectroscopic signals revealing the SN → TT1 SF in real time. We demonstrate that SF can proceed directly from the upper state SN, bypassing the lowest excited state, S1. We determine the main SN → TT1 reaction pathways and show by computer simulation and spectroscopic measurements that the SN-initiated SF can be faster and more efficient than the traditionally studied S1 → TT1 SF. We claim that the SN → TT1 SF offers novel promising opportunities for engineering SF systems and enhancing SF yields.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals.

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Engineering Cavity Singlet Fission in Rubrene.

By employing the numerically exact multiple Davydov D2 ansatz, we study cavity-manipulated singlet fission that is mediated by polaritonic conical intersections for both one- and two-molecule systems. The population evolution of the TT state and the cavity photons is carefully examined in search for a high fission efficiency via cavity engineering. Several interesting mechanisms have been uncovered, such as photon-assisted singlet fission, system localization via a displaced photon state, and collective enhancement of the fission efficiency for the two-molecule system. It is also found that the system localization process in the two-molecule system differs substantially from that in the one-molecule system because of the appearance of a novel central polaritonic conical intersection in the two-molecule system. It has been demonstrated that the cavity-controlled singlet fission process can be switched on and off by controlling the average pumping photon number.

Accurate Simulation of Spectroscopic Signatures of Cavity-Assisted, Conical-Intersection-Controlled Singlet Fission Processes.

A numerically accurate, fully quantum methodology has been developed for the simulation of the dynamics and nonlinear spectroscopic signals of cavity-assisted, conical-intersection-controlled singlet fission systems. The methodology is capable of handling several molecular systems strongly coupled to the photonic mode of the cavity and treats the intrinsic conical intersection and cavity-induced polaritonic conical intersections in a numerically exact manner. Contributions of higher-lying molecular electronic states are accounted for comprehensively. The intriguing process of cavity-modified fission dynamics, including all of its electronic, vibrational, and photonic degrees of freedom, together with its two-dimensional spectroscopic manifestation, is simulated for two rubrene dimers strongly coupled to the cavity mode.

Dynamics of disordered Tavis-Cummings and Holstein-Tavis-Cummings models.

By employing the time-dependent variational principle and the versatile multi-D2 Davydov trial states, in combination with the Green's function method, we study the dynamics of the Tavis-Cummings model and the Holstein-Tavis-Cummings model in the presence of diagonal disorder and cavity-qubit coupling disorder. For the Tavis-Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green's function method to corroborate numerically exact results of Davydov's Ansätze. For the Holstein-Tavis-Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.